[LONGITUDINAL]
input = {Cc, Emax, ED50, gamma, E0}
Cc = {use=regressor}

EQUATION:
Cc_sat = max(Cc, 0)
if Cc_sat <= ED50
  r = (Cc_sat/ED50)^gamma
  cc = r/(1+r)
else
  r = (ED50/Cc_sat)^gamma
  cc = 1/(1+r) 
end
A = Emax * cc
S = E0
E = A + S

OUTPUT:
output = E